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1 Introduction 

The multidimensional gravitation theories are very attractive in the context of the unification of 
fundamental interactions. Moreover, several modern theories require space-time to have more than 
four dimensions [1-8]. The nonobservability of additional dimensions in such theories needs an 
explanation. Among different possible ways of such explanation the hypothesis about dynamical 
contraction of internal manifold during expansion of the universe is very popular. This idea is 
realized in many exact cosmological solution of multidimensional Einstein's equations [9-22]. As a 
rule such models require additional fields and do not avoid initial big bang singularity. The intro- 
duction of additional fields in multidimensional gravitation theories destroy their pure geometrical 
character and require an additional motivation [4] . Such motivation may be done in the framework 
of some generalizations of Riemannian geometry. In four dimensional case such generalization in 
several cases leads to removing of cosmological big bang singularity [23-25] . That is why the unifi- 
cation of generalized geometric structures and multidimensional gravity seems to be very attractive. 
Unfortunately, only in several papers the multidimensional gravitation theory and cosmology are 
considered in the scope of some generalization of Riemannian geometry [26-28]. 

One of the simplest generalization of the Riemannian geometry is the integrable Weyl geometry 
with the connection components 

where T^^ are the Christoffel symbols, uja = a;,^, a; is a scalar field, 5p are the Kroneker symbols, 
ga/3 is a metric tensor; the small Greek indices take values from to n — 1, n is a dimension of 
space-time. The Ricci tensor and the curvature scalar of the connection (1) are equal to 

„ ~n — 2 l_n — 2/ x\ 

R/xu = Rixu H l^^i^Wv + 2^/^'^^ ~ y^i^'^v - 9iiv^ ^\ ) , (2) 

i? = ^+(n-l)Da;-^^?— ^fc^wV, (3) 

where the tildes denote the quantities calculated in the connection F^^ , two parallel vertical bars and 
□ denote the covariant derivative and the d Alembert operator of this connection. It is necessary to 
note that the integrable Weyl space-time is also conformally-Riemannian, since there is a conformal 
transformation of metric tensor g^p which maps the Riemannian space-time into integrable Weyl 
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space-time. As the integrable Weyl space-time is defined by the pair {ga/s, the gravitation theory 

in this space-time docs not coincide with Einsteinian general relativity because the field uj must 
be contained in the Lagrangian independently from g^p and cannot be excluded by the conformal 
transformation. 

Some features of the Einsteinian cosmological models with scalar fields were recently consid- 
ered by several authors [9,11,19,21,28,29-37] both in 4-dimensional and in (4-|-d)-dimensional space- 
times. The cosmological models in four-dimensional Wcyl-integrable space-time were recently con- 
sidered by Novello et al. in [23], where the existence of nonsingular open cosmological models was 
shown. The appearance of Weyl geometry in multidimensional cosmology was discussed also in 
[27]. 

In this paper we consider the influence of Weyl geometry on the evolution of Friedman- Robertson- 
Walker (FRW) cosmological models in multidimensional gravitation theory. As usually the space- 
time is assumed to have the structure of direct product x of four-dimensional FRW space-time 
and d-dimensional interior space that is supposed to be d-sphere S"^ or c/-torus T*^. The 
metric of space-time is supposed to be block-diagonal 

ds"^ = dt^ - a^{t) i^J^^ + r^c^O^j - g^i,du''du\ (4) 

where k = +1, 0, —1 for closed, plane and open models, dO,"^ is a line element on two-sphere, n", 
a = l,...,d, and widega/3 are the coordinates and metric tensor of the interior space V^. Once 
we consider only spatially homogeneous FRW cosmologies, it is natural to make the Weyl scalar 
field u; to be a function of cosmic time t only: lo = Lo{t). We consider both vacuum case and 
non vacuum case with the additional scalar field (p with non minimal coupling. The 4-dimensional 
case will be briefly considered also for completeness. The existence of the conformal map between 
Riemannian and integrable Weyl space-times may be used for generation of exact solutions from 
the known solutions of general relativity. Such approach admits obtaining only the particular 
solutions. Therefore to demonstrate general qualitative behavior of the models we solve the system 
of cosmological equations numerically with initial values given at t = and satisfying the constraint 
equation. For that purpose we use adaptive numerical methods with automatic choice of integration 
step and with the stiffness checking. The geometrical units where G = c = 1 are used in what 
follows. 



2 Integrable Weyl cosmology in vacuum 

Following [23] we shell consider the vacuum cosmological models in the gravitation theory with the 
Lagrangian 

L = R + ^LOaUj'^ (5) 

where R is defined by (3) and ^ = const. After excluding the total derivatives of the scalar field 
Lagrangian (5) takes the form 

L = il-'"-'><7^'-'^.»., (6) 

So, the theory differs from the Einstein theory with the massless scalar field by the coefficient 

before the square of the scalar field gradient and has different geodesic lines. Note also that due to 
the definition of the Weyl connection (1) the scalar field uj cannot be renormalized and hence the 
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coefficient ^ before w'^Wq cannot be put to ±1 as it may be done in the pure Einstein theory with 
massless scalar field. Variation of (6) with respect to the pair (^0/3, a;) of independent variables 
yields the equations 

~ 1 ~ (n-l)(ra-2)-4^ / 1 „ \ . 

Rixv - l^QtiuR 1 U^iii^y - ^gfj.iy(^ I = 0, (7) 

and 

□cj = (8) 

The equations (7), (8) coincide with the Einstein equations for the massless scalar field, whose 
solutions for the FRW cosmological models were investigated both in four-dimensional [32] and 
multidimensional cases [37]. By this reason here we only summarize briefly the main results. 

2.1 Four- dimensional case. As the scalar field w is a function on t only, equation (8) yields the 
first integral 

(9) 

where ovcrdot denotes time differentiation and 7 = const is the integration constant. Due to (9) 
equations (7) take the form 

a^ + k-^ = (10) 

and 

2ad + a^ + k + -rhr = (11) 

where A = (3 — 2^). As it is easy to see from (10), only singular and static solution of equations 
(lO)-(ll) exist if A > 0. For negative values of A solution exists only for the open models. In 
this case a{t) > qq = ~ 3)7^/12 and so the cosmological singularity is absent. The qualitative 
behavior of scale factor a{t) for negative A is shown at figure 1 and its features are discussed in 
detail in [23]. 

Fig. 1. 

Qualitative behavior of the scale factor a{t) of the open universe in four-dimensional 

Weyl-integrable space-time. 

2.2. Multidimensional case. In the multidimensional case the behavior of the model depends not 
only on the parameter ^, as in the previous case, but on the structure of the interior space also. For 
simplicity only 5- and 6-dimensional models will be considered in the following. We consider these 
two cases separately. The main qualitative features of models in general n-dimensional (n > 6) 
case are the same as in 5- and 6-dimcnsions. 

2.2.1. 5-dimensional models. In 5-dimensions space-time interval (4) reads 

/ dr^ \ 
ds^ = dt^ - a^{t) YZTj^ + ^^^^"^ - s^{t)du^ (12) 

where u G 5^ is the interior space coordinate. Assuming, as above, a scalar field a; to be a function 
of the cosmological time only, the first integral of equation (8) takes the form 
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where 71 = const. Due to (12)-(13) equations (7) become after simplification 



a \a/ as a 

and 



ds /d\2 3k 7?(3-£) ^ 
3— + 3 - + -iii^=0, 14 
as \aj za^s^^ 

a f a s k ^ 

7 + 2(7) +77 + 72 =0' (15) 



The last equation has the first integral 



s ^a s ^ 

- + 3 — = (16) 
s as 



»=3 (iv) 



where 72 = const. It is easy to see that analogous to the four-dimensional case the nonsingular 
solutions of equations (15), (17) exist only for the open models [k = —1). In this case for i < the 
scale factor of 3-spacc a(t) decreases monotonically from infinity to its minimal value qq and then 
grows to infinity at t > 0, while the Weyl field uj{t) and the scale factor of interior space evaluates 
monotonically from U- = limj_>_oo w(t) and s_ = limj_>_oo s(i) to s+ = limj_>oo s(t), where oq, 
L0± and s± are defined by the integration constants and may have arbitrary values. Note that if 
72 < than the constants s_ and s+ satisfy the condition s_ > s+ and so the standard dimensional 
reduction scenario is realized. The typical shape of the functions a{t) and s{t) are shown on the 
figures (2.a,b). 

Fig. 2. 

The qualitative behavior of the scale factors a{t) (2a) and s{t) (2b) of five-dimensional Weyl 
integrable cosmological model with open 3-space. 

The figure (2a) shows that unlike the 4-dimensional case the evolution of 3-space in 5-dimensional 
model is time- asymmetric. This asymmetry appears because the equation (15) depends not only 
on s{t) but also on the time-asymmetric interior space scale factor s{t). 

2.2.2. 6- dimensional models. In 6-dimensional case we consider two types of topological struc- 
tures for the interior space: the 2-sphere and 2-dimensional torus T^. Therefore the space-time 
metric (4) may have one of two forms 

+rW]-s\t)(^^ + uV] (18) 




or 

, ^ + r'^dn'^ - sj(t)du^ - sl(t)dv'^, (19) 
1 — fcr^ J 

where {u, v} are the coordinates on S'^ or respectively. First integrals of equation (8) take the 
form 



for metric (18) and 
for metric (18). 



^ = ^^. (21) 
a'^sisz 
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Equations (7) for the metric (18) after simplification take the form 



+ 2(t) +277 + -2=0' (22) 



a \a/ as 



s /s\^ „ds 1 



-+ - +3— + ^ = 0, (23) 



and the constraint equation 



„a / a „s\ 3k 1 oi (5 — £) , 
3- - + 2- + — + — - ^ „ = 0. 24 

The first two equations are dynamical and the last is the constraint. 

It is easy to see that only singular solutions of equations (22)-(24) exist: the scale factor s(t) 
of the interior space evolves from zero at t = to its maximal value Smax and return to zero at 
t = t\> to. The behavior of a(t) depends on the sign of k. Namely, if fc = +1 then the qualitative 
evolution of a{t) is the same as the evolution of s{t). If A; = than a{t) increase from zero at t = to 
to infinity at t = ti or decrease from infinity to zero; the unstable solutions with a{t) = cOnst are 
also exist. Finally, if /c = — 1 then a{t) evolves from infinity at t = to to its minimum amin and 
then grows to infinity at t = ti. 

For the metric (19) equations (7) after simplification read 

a a / si S2\ „/a\^ 2fe „ ,„_, 
- + - - + - +2 - +^ = 0, 25 
a a \si Si) \a J a 

51 0, Sl Sl So ^ , . 

— + 3 — - + —— = 0, (26) 
si a si si S2 

So d So Si So , , 

— + 3 — ^ + — — = 0, (27) 

52 a S2 Si S2 

and the constraint equation 

3^f^ + ^ + ^)+^^ + ^-4^x# = 0- (28) 
a \a si 52/ si S2 za^spg 

As in 4- and 5-dimensionaI cases the nonsingular solutions of the equations (25)-(28) exist only 
for the open models (k = —1). Analogously to 5-dimensional case the scale factor of 3-space a{t) 
in these models decreases monotonously from infinity to its minimal value uq and then grows to 
infinity at t — +oo, while the scale factors Si(t), z = 1, 2, of interior space changes monotonously 
from Si- = lim(^_(-o Si(t) to Sj+ = limt_>oo Si(t). The necessary condition for the realization of the 
dimensional reduction scenario in this case are defined by the following inequalities 



9i(5-0 

% 

and 



- + > 0> (29) 

05 2a3sioS2o 



si(0) <0,S2(0) <0 (30) 

It is necessary to note that inequality (29) is the necessary condition for si and S2 to be of the 
same sign. The time behavior of scale factors a{t), Si{t) and S2{t) in this case is qualitatively the 
same as in 5-dimensional case (Figure 2). 
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3 Integrable Weyl cosmology in theory with non minimal scalar 
field 



In this section we consider cosmological models in gravitation theories with Lagrangian 

L = r(i+ . ^ ^y )+^Lo''LOa + ri<p''ipa, (31) 



2(n - 1) 

where R is defined by (3), 99 is a real scalar field, r] = ±1 and = const as above. In the limiting 
case if = const Lagrangian (31) coincides with (5) while in another limiting case lo = const it 
coincides with the Lagrangian for the conformal-invariant scalar field. 
The substitution of (3) into (31) gives after simplification 

L = R^1 + - w'-^e. - ^"""^^^Y^^"^^ ^"^" " ^^^V^u-u^ + (32) 

where the total derivatives of the scalar fields are omitted. 

Variation of (32) with respect to independent variables g^i,, to and cp yields the equations 

f; ^ ?C\f. \ (n- l)(n-2) -4^ / 1 
R^^u - -^Qi^uR J 1 1 + 2{n-l) ) 4 V"^'' ~ 2^^"^ 

1, . n — 2 2 f ^ n \ ^ f ^ \ 

(^v - -^^^ + g^^i^'p'^Va ~ 2 ) " ^' ^^^^ 

/(n- l)(n-2) -4^ n-2 o\„ „ „ 
'-^ ^ + □a;-^D(^-(^„(^,""=0, (34) 

and 

/ 1 ~ n — 2 \ 

r)D^p -{Dlo + -R —uj," u;,y]ip = 0, (35) 

V n — 1 4 / 

Equation (35) shows that non-Riemannian nature of space-time geometry in the considered 
model leads to the effective mass generation for the scalar field ip. 

3.1. Four- dimensional models. In four-dimensional case the equations (33)-(35) consist of the 
constraint equation 



-^^-2 1 ■■ 1 2-2 3-2^. 



and three dynamical equations 



a 1 .. faV I (p'^\ 5d . 1 .0 (p'^\ 2k 



and 



\^^-2^ + 3jQ-ipyi+(9-6^ + ^(^^^ ^d; - 3%<p - = 0, (38) 

d o. . ^ 1 ■ 2 ( 2k 

■qip — ipu + 2ip- + 3- [nip — (pu) + -ipu +2(p { - ] + —^ip = 0. (39) 
a a 2 \aj a"^ 
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The coefficients before a/a, oj and (p in the equations (37)- (39) depend both on the parameters 
^, rj and on the scalar field ip. The determinant of the matrix of coefficients before a/a, u and (p is 
equal to 

d = (2 - I) (^4 + (277 + f - M _ 4) (^2 + 6r/ - 47?e 
The points where d = are the singular points of the system (37)-(39). These points are not 

described by the system (36)- (39) because for fixed rj and ^ equation d = const defines not more 

than four fixed values of ip and the system (36)- (39) reduces to the first order system. Therefore 

the initial value of the field (f must be from the open set d 7^ 0. 

For 77 equation d = divide the half-plane (^, y?^), on three regions that will be denoted as A, 

B and C, while for 77 = —1 there are only two regions A and B (figure 3a,b). The behavior of the 

model depends on the region where the point (^, (Pq) is situated. 

Fig. 3. 

The set d = for 77 = 1 (3a) and 77 = —1 (3b) in 4-dimensions. 

Numerical investigation of equations (31)- (33) shows that for the closed {k = 1) and flat {k = 0) 
cosmological models only singular solutions exist for any initial conditions. For the open models 
(fc = —1) if the pair (^,1^0) defines the point in the region B (both for 77 = 1 and 77 = —1) or 
C (for 77 = 1) than only singular solutions of the equations (37)-(39) exist. If the pair {C'Po) 
defines the point in the region A then solutions may be both regular and singular. The numerical 
investigation does not permit to find the exact conditions of regularity, but it shows that both 
regular and singular solutions are stable against finite perturbations of the initial conditions. The 
typical qualitative behavior of the universe scale factor a{t), Weyl field ui and the matter scalar 
field if are shown at figure 4a-c. 

Fig. 4. 

The qualitative behavior of the universe scale factor a{t) (4a), Weyl field lo (4b) and the matter 
scalar field (p (4c) in the open four-dimensional Weyl-integrable space-time model. 

The universe scale factor a(t) in the typical nonsingular solution evolves from infinity at i = — 00 
to its minimal value oq = a(0) and then grows to infinity at t — > 00 (figure 4a). Both scalar fields, 
the Weyl field uj and the field (p evolves between two limiting values: from uj- = ]imt^-^uj{t) and 
</?_ = limt_»_oo to uj^ = limt^^oo i^(i) and 99+ = liuit^oo The difference in the evolution of 

these fields is that the field lo evolves monotonously (figure 4b) while the field ip near t = (i. e. near 
the minimum of a{t)) may have several intermediate extrema with one absolute maximum if r/ = 1 
(figure 4c) or absolute minimum if 77 = —1. As p{t) for big \t\ tends asymptotically to constants, the 
model evolves asymptotically as an empty Weyl cosmological model that is considered in section 
2.1. It is necessary to note also that the evolution of the universe scale factor a(t) has a small 
time-asymmetry in comparison with the case of the empty space. This asymmetry is a result of 
non symmetrical evolution of the matter field (p because the field equations (37)- (39) contain both 
(p and if). 

3.2. 5-dimensional models. In 5-dimensional case equations (33)-(35) after simplification be- 
come 
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J + 2- + Ui + ^l + ^^ =0, 41 

{Lp^\ [a a f a s\ 2k^ ioip ii>(p fas (3 \ , , 

-(^-2« + 6)-w + -(3^ + i) (3^-25 + 6) -w(3^ + i + ^)=0. (43) 

and 

(/? /3d s\ .. .. . / 3lL' d s\ . f „a s\ 3ix>a fa s\ 3kip 

The equation (40) is the constraint that must be satisfied by the initial conditions and the 
equations (41)-(44) are the dynamical. The determinant of the matrix of the coefficients before 

d/a, s/s, Co and (p in the dynamical equations (41)-(44) is equal to 

The qualitative features of function d(^, r/, ip) are the same as in 4-dimensional case: for r/ = 1 
equation d = divides the half-plane (p^ > 0) on three regions that are denoted as A, B and 
C while for r/ = — 1 there are only two regions A and B (figure 5a,b). The behavior of the model 
depends on the region where the point (^, (/Pq) is situated. 

Fig. 5. 

The sets c? = for = 1 (5a) and rj = —I (5b) in 5-dimensions. 

Numerical investigation of equations (40)-(44) shows that as well as in the previous 4-dimcnsional 
case only singular solutions exist at any initial conditions for the closed {k = 1) and flat (k = 0) 
cosmological models. For the open models (fc = —1) if the pair {(,,Pq) defines the point in the 
region B (both for 77 = 1 and 77 = — 1) or C (for rj = 1) than only singular solutions of the equations 
(40)-(44) exist, while if the pair (Cpq) defines the point in the region A than the solution may 
be both regular and singular. The regularity of solutions depends on the constants of integration 
that may be considered as the initial conditions at t = 0. It was found that the regularity of 
solutions depends mainly on the signs of s(0), u;(0) and i^(0). Their possible combinations that 
give nonsingular solutions of equations (41)-(44) are represented in table 1. The last column of 
this table shows the general direction of the interior space evolution by means of the signs of the 
difference A = s+ — s_, where s± = lim.t^±oo s{t). 

Table 1. 

Conditions of the solutions regularity and the direction of s{t) evolution 
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signs (0) 


sienwfO) 


signwfO) 




-1 


-1 





-1 


-1 


+1 





-1 


-1 


-1 


+1 


-1 


-1 


+1 


-1 


-1 





+1 





-1 





-1 





+1 


+1 


+1 


-1 


+1 


+1 


-1 





+1 


+1 


+1 





+1 


+1 


-1 


+1 


+1 



The typical behavior of the nonsingular solution of the equations (41)-(44) for 77 = 1 is shown 
at the figures 6a-d for the case A < 0, i. e. for the contracting interior space. 

Fig. 6. 

The qualitative behavior of the scale factors a{t), (4a) and s(t) (4b), the matter scalar field ip{t) 
(4c) and Weyl field u!{t) (6d) in the nonsingular 5-dimensional Weyl-integrable space-time model. 

In general nonsingular solution the radius of the universe changes monotonously from infinity at 
f = — 00 to minimal value oq and then grows to infinity (figure 6a), while the radius of the internal 
space starts from s_ = limt_»_oo s(t), passes through several (one or two) intermediate extrema, 
that are situated near minimum of a{t) and may be absent in some cases, and then changes to 
s+ = lim^^oo •s(^) (figure 6b). Note that s_|_ and s_ may be of the same or different order. The 
field (p evolves analogously to 4-dimensional case (figure 6c). Note that the extremal points of the 
functions a{t), s{t) and ip{t) do not coincide with each other in general case and the function a{t) 
is time asymmetrical especially near its minimum. Finally the Weyl field uj changes monotonously 
between two limiting values: a;_ = lim^^-oo Ci;(t) and a;+ = lim.t^oo'^i.t) (figure 6d). In the case 
rj = —1 the model evolves as above but the extremal points of the field change type: minimum 
become maximum and vice versa. 

4 Concluding remarks 

We have considered the qualitative evolution of multidimensional cosmological models based on the 
integrable Weyl geometry both in vacuum space-time and in the presence of nonminimal scalar field. 
The existence of nonsingular solutions of field equations for open cosmological models that realized 
the dimensional reduction scenario was demonstrated. It was shown that in multidimensional case 
the evolution of the scale factor of the universe a{t) becomes time-asymmetric unlike the four- 
dimensional case. We have shown also that all nonsingular cosmological models considered above 
have some common features. In particular the evolution of the universe scale factor (radius) a{t) 
for big \t\ is asymptotically linear. Further in all nonsingular models Weyl scalar field ui{t) as well 
as the matter field (p{t) in the models with nonminimal coupling tend asymptotically to constants. 
So the models tend to the pure Einsteinian models of the corresponding dimensions and the change 
of the collapse era into expansion one may be considered as a cosmological phase transition induced 
by the transition of scalar fields io(t) and ip{t) from one stationary state uo = uj- and tp = ^p- into 
another stationary state a; = a;+ and ip = ip^. At the late stages of the universe evolution the fields 
a;(t) and <^(t) are unobservable. 
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There are several qualitative differences between the vacuum models and the models with non- 
minimal scalar field. First of all in vacuum models the existence of cosmological singularity depends 
only on the parameters of the theory while in the case of nonminimal scalar field it depends on the 
initial conditions also. Secondly, in the models with nonminimal scalar field the evolution of the 
internal space scale factor s{t) may be nonmonotonous. In the typical scenario one of the limiting 
values of s{t) at t = ±00 is much smaller than another but in several models both limiting values of 
internal radius ,s{t) may be arbitrary small and it become finite only near minimum of the universe 
scale factor a{t). 

We have discussed here only the models with the one- or two- dimensional interior space because 
if interior space has dimension d > 3 and direct product topology of torus on several spheres then 
the models have the same qualitative features as considered above. In particcular, the nonsingular 
solutions exist only for toroidal interior space topology. 

The models considered above show that the real geometrical structure of space-time may have a 
non-Riemaniann nature but the universe may evolve in such a way that its non-Ricmaniann nature 
is essential only near t = and become unobscrvable at late stages of the evolution. Therefore, 
the consideration of generalized geometrical structures in multidimensional cosmology may be of 
a considerable interest. In particular, the models considered above may be generalized in the 
following manner. First of all, both Weyl scalar field uj(t) and matter field (p{t) may be massive 
and have nonlinear potential. Secondly, the possible influence of the cosmological term A must be 
considered also. At last, the term Rip'^/2{n — 1) in the action integral (25) may have negative sign. 
One may suppose that in this case nonsingular solutions of the field equations may be obtained 
not only for open models, but for closed and flat models also. These possibilities will be considered 
elsewhere. 
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